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Abstract 

We compute finite- volume corrections to nucleon matrix elements of the axial- vector current. We 
show that knowledge of this finite-volume dependence — as well as that of the nucleon mass — 
obtained using lattice QCD will allow a clean determination of the chiral-limit values of the nucleon 
and A-resonance axial-vector couplings. 



1 



I. INTRODUCTION 



The nucleoli axial charge, qa, is a fundamental quantity in QCD as it in some sense quantifies 
spontaneous chiral symmetry breaking in the low-energy hadronic theory. It is known very 
accurately from neutron beta decay experiments, = 1.2670 ± 0.0030 (in units of the 
vector charge gy) 0] and therefore serves as an important test of nonperturbative, first- 
principles calculations of hadronic properties using lattice QCD. Fortunately, it is relatively 
straightforward to compute gA in numerical lattice QCD simulations. In spite of this, there 
is still no consensus as regards 04 from lattice QCD 0, S 0, S 0, 0, H 01- In fact, a trend 
toward under predicting [2, lEJi h as led to some controversy regarding the possibility 
of large finite- volume effects |lOl A recent quenched simulation using domain- wall 

fermions over several volumes finds that large finite-volume effects do seem to account for 
the discrepancy The somewhat tentative current state-of-affairs is primarily due to the 
fact that current computational limitations use quark masses, m^, that are significantly 
larger than the physical values, lattice spacings, a, that are not significantly smaller than 
the physical scales of interest, and lattice sizes, L, that are not significantly larger than the 
pion Compton wavelength |12|. It is confidence in the extrapolations of these quantities 
that will allow a confrontation between lattice QCD predictions for gA and other hadronic 
observables and experiment. While the dependence of gA on the lattice parameters can be 
described by an effective field theory (EFT), calculability requires maintaining the hierarchy 
of mass scales: |p| , m.,^ -C <^ a^^, where |p| is a typical momentum in the system 
of interest, is the pion mass and A^ ~ 2\/2nf is the scale of chiral symmetry breaking 
(/ = 132 MeV is the pion decay constant). Lattice simulations are only now beginning to 
achieve the hierarchy of scales necessary to utilize a perturbative extrapolation. 

Here we will be concerned primarily with the finite-volume dependence of gA- In a 
spatial box of size L, momenta are quantized such that p = 27rn/L with n G Z. The 
EFT momentum hierarchy then requires maintenance of the additional inequality f'L 3> 
1. This bound ensures that (non-pionic) hadronic physics is completely contained inside 
the lattice volume. In addition, the bound {mT^LY{fVf' 3> 1 ensures that the lattice 
volume has no effect on spontaneous chiral symmetry breaking [isL l3 | . These two bounds. 



taken together, then imply that in order to have a perturbative EFT description mT^L^ 1. 
When {ttIt^LY {f L^S. 1, and therefore m^L ^ 1, momentum zero-modes must be treated 
nonperturbatively 0, 3| and one is in the so-called e-regime ^. 



We will consider the range of pion masses ^, 130 MeV< 300 MeV, and therefore we 
will take L> 2 fm, keeping in mind that the EFT may be reaching the limits of its validity 
when this bound on L is saturated, particularly when the pions are light. For the observables 
considered here, finite-volume effects tend to be small for L > 4 fm. It is therefore of 
interest to have control over the finite-volume dependence of hadronic observables in the 
range 2 fm < L < 4 fm. Chiral perturbation theory (^PT), which provides a systematic 
description of low-energy QCD near the chiral limit, is the appropriate EFT to exploit the 
hierarchy of scales described above and to describe the dependence of hadronic observables 
on L llil. ITgI 17. l^L Recent work has investigated the finite- volume dependence in the 



meson 



cM. iia 11 g liQll. Kecent worK nas mvestigatea tne nnite-voiume depenaence m tne 
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The chiral-limit considerations of Ref. fall in the e-regime. However, to our knowledge no system- 

atic finite-volume calculation of baryon properties has been done in the e-regime. 
^ The current upper limit of this range has been estimated recently by one of the authors [l5j 
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sector. 

In this paper we compute the leading finite- volume dependence of the axial- vector charge 
of the nucleon in heavy-baryon xPT (HBxPT), including the A- resonance as an explicit 
degree of freedom. The finite-volume corrections to the axial-vector charge of the nucleon 
depend on the A-nucleon mass splitting and on the chiral-limit values of the nucleon, A- 
nucleon and A axial-vector charges. Traditionally, the nucleon and A axial couplings have 
been estimated using the spin-fiavor SU (4) symmetry of the quark model, and in recent 
work [3^ the authors have conjectured the chiral-limit values of these couplings. We point 
out that lattice QCD measurements of finite-volume effects in the axial-vector charge (and 
mass) of the nucleon will provide a clean determination of the nucleon and A-resonance 
axial-vector couplings. 



II. THE NUCLEON AXIAL CHARGE IN A FINITE VOLUME 



At one-loop level, the matrix elements of the axial-vector current between nucleons of 
flavor "a" and "b" may be written as 



{Nb\ J^,5 \Na) = [ Tab + Cab ] 2UbS^Ua , 



(1) 



where Cab represents a counterterm with a single insertion of the light-quark mass matrix. 
The leading-order Lagrange density describing the interactions between the pions and the 
low-lying baryons is 



C = 2gA NS^A.N + g^N 



-pr^abcM Arl T^T 1 

T A'^auNbecd + h.c. 



+ 2gAA T^S^^A^n 



(2) 



This Lagrange density gives rise to the diagrams in Fig. ^ which are the leadin g on e-loop 
contributions to the axial-current matrix elements. In the isospin limit one flnds [35| 



NN 
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451^ Ji(m^,0,/i) + 4 ( ^(ijvfl'A + — S'ATvfi'AA ) Ji(m^,A,/i) 



3 32 
+ -^9a -9\n9a A^i(m^, A,/i) 



(3) 



where Ji(m, A, /i), Ri{m, /i) and Ni{m, A, fi) are loop integrals defined in the Appendix and 
A is the A-nucleon mass splitting. All r[e] poles have been subtracted. They — and their 
associated counterterm c^n — need not concern us here as the finite- volume corrections do 
not depend on the ultraviolet behavior of the theory at leading one-loop order. All of the 
couplings (including /) in eq. take their chiral-limit values. 

Using the notation 6l = "^iL) — ^{oo) to denote the finite-volume corrections to the 
quantity and the results obtained in the Appendix, the finite- volume corrections to Fattv 



are 



5l (Tnn) = SgA 



mz 



25 

^^Fi + ( 9'an9a + —9an9aa ) F2 + ^^Fg + g'^^gAF^ 



,(4) 



where 



Fi(m,L) = Y: 



-K'o(mL|n| 



Ki{mL\n\) 
mLlnl 
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FIG. 1: One-loop graphs that contribute to the matrix elements of the axial current in the nucleon. 
A solid, thick-solid and dashed line denote a nucleon, a A-resonance, and a pion, respectively. The 
solid-squares denote an axial coupling given in eq. while the crossed circle denotes an insertion 
of the axial-vector current operator. Diagrams (a) to (e) are vertex corrections, while diagrams (f) 
and (g) give rise to wavefunction renormalization. 
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,(5) 



and Ka{z) is a modified Bessel function of the second kind. Tlie extension of this result to 
PQQCD, including strong isospin violation, is straightforward to extract from Ref. [s^ using 
the results of this paper. We do not give an asymptotic expression for 6gA as we do not find 
it useful for L < 10 fm for the pion masses of interest, however, it may be found by taking 
the appropriate asymptotic limits of eq. (0) using technology developed in Ref. [I^li^. One 
sees in eq. © that, as itIt^L 0, the Fi diverge, signaling the transition to the e-regime 
and the necessity of a non-perturbative resummation. 
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III. EXTRACTING AXIAL CHARGE FROM LATTICE QCD 



The finite- volume corrections toT]\iN depend only on infrared quantities, i.e. the axial- vector 
charges and the pion decay constant, the meson mass and the A-nucleon mass-splitting. 
Hence, with precise determinations of / (chiral-limit value), and A, lattice data at 
several different values of L will allow a determination of the axial- vector charges. However, 
in order to separate the various contributions to eq. (@)), one must ensure that the Fi scale 
differently over the relevant values of L. In Fig. El we plot Fi and the ratios F2/F1, F3/F1 
and F4/F1 as functions of L for various pion masses. For F2 and F4 we use A = 293 MeV. 
It is clear from Fig. |21that the ratios of the Fi scale differently and therefore, in principle, 



0.3 




L (fm) L (fm) 



FIG. 2: Plot 0/ Fi and the ratios F2/F1, F3/F1 and F4/F1 vs. L. The solid and dashed lines 
correspond to = 139 MeV and 300 MeV, respectively, for the physical A-nucleon mass splitting, 
A = 293 MeV. 



the coefficients of the Fi in eq. (jH) r nay be extracted from the L dependence of Sqa - 

In a recent paper by the authors [3J|, based on earlier work by Weinberg 37, 3^ 3^, 4^ . 
it was conjectured that in the chiral limit, the helicity one-half components of the nu- 
cleon, A and the Roper (A^'(1440)) fall into the reducible (2,3) © (1,2) representation of 
SU (2) L SU {2)fi with maximal mixing. Denoting the mixing angle between the irreducible 
representation as ip (with maximal mixing corresponding to ip = vr/4), the conjecture deter- 
mines the chiral-limit values qa = i + (2/3) cos^ ipy 9 an = —2 cosip and (^aa = —3. Inserting 
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these values into eq. (jH) leads to 



5 9/ 



m: 



Fi 



2Fi + — F2 + ^F3 + 4F4) cos^V^ 



+ - (F1 + 2F2 + 2F4) cosV + ^Fi cos^V^ 



(6) 



It would be interesting to have a direct lattice determination of ip using this formula. The 
spin-flavor SU{4) (naive constituent quark-model) results are recovered with ip = 0. How- 
ever, the conjectured values (with ip = tc/A) are in much better agreement with existing 



experimental knowledge [3^ |41[ . We use eq. © to estimate our current knowledge of the 
finite- volume dependence of the nucleon axial- vector charge. This expression is plotted as a 
function of L for various pion masses in Fig. 01 for the two cases ip = ir/A and ip = 0. Varia- 
tion of ip provides a measure of the experimental uncertainty associated with the chiral-limit 
values of the axial- vector couplings ^. It is encouraging that the two scenarios lead to quite 
distinct predictions for Sqa, and therefore a precise determination of the volume dependence 
of gA will allow for a determination of the mixing-angle ip. In both cases it is clear that for 
L> 2 fm, finite-volume effects are at the few-percent level for all relevant pion masses. 
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FIG. 3: The volume dependence of qa for chiral-multiplet mixing-angles ip = 7r/4 and ^ = 0. The 
left panel shows 6gA vs. L with ip = 7r/4, where the solid, dotted and dashed lines correspond 
to = 139 MeV, 200 MeV and 300 MeV, respectively. The right panel shows 5gA vs. L with 
Ip = (spin-flavor S'[/(4) values of axial-vector couplings). The physical A-nucleon mass splitting, 
A = 293 MeV, is used for both panels. 



IV. CONCLUSIONS 

It has long been known that the infinite-volume S-matrix can be extracted from power- 
law suppressed finite- volume effects that arise when a two-particle state is put in a finite 



^ For a recent discussion of current knowledge of the chiral-limit value of gA, see Ref. [l 
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volume ji^, 3, 3| ) cind very recently it has been shown that this method may be ex- 



tended to include the effect of external electroweak gauge fields j43[. Therefore, a lattice cal- 
culation of the energy-levels of a pion and a nucleon in a finite volume can, in principle, allow 
for an extraction of the axial-vector couplings. Important information may also be extracted 
from exponentially-suppressed finite- volume effects that arise from quantum-loops 3|- An 



important observation is that ultraviolet physics (counterterms) enters the chiral expansion 
for finite- volume effects beyond leading one-loop order in the expansion • Finite- volume 
effects therefore offer a clean probe of infrared physics. Moreover, this method is optimal 
for ttItt < A where the A-resonance is unstable and a direct probe of A properties is prob- 
lematic. To conclude, we have computed the leading finite-volume corrections to nucleon 
matrix elements of the axial-vector current and shown that a lattice QCD measurement 
of this finite-volume dependence can determine the chiral-limit values of the axial-vector 
charges of the nucleon and A-resonance. 
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APPENDIX: LOOP INTEGRALS 

In this appendix we review some standard one-loop integrals that arise in HBxPT 1^ and 
give their finite-volume dependence. First we consider the generic one-loop integral 

1 1 



/o(m,A,/i) = At^y 



(27r)" qo — A + ie — rn? + ie 



87r2 

i 



A log — - 2 A - VA2 -m'^ + ie log , , 

vr^(m,A,/i) , (A.l) 

where 7rJ-'{m, 0, fi) = Tcm, e = 4—n and we have subtracted the Ar(e) divergence. Evaluating 
the energy integral yields 

i f°° P f d"'~^q 1 

JoK A,/.) = d(3 J ]3/. ' (A.2) 

where /9(A)^ = + 2AA + and we have performed a change of variable that is valid only 
for A > 0, as the relation is non-invertible for A < 0. Using the master relation 



2 



'A - ^/A? ~m? + ie 



which has been derived previously j30|, |33[ , one finds the finite- volume corrections 

i roo Q 

(/„(™, A. ,)) = ^ £ d3 E K,(3L\n\) • (A.4) 
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Notice that there is no renormahzation-scale dependence. In general, the integral over (3 
cannot be performed analytically, however, for A = 



5l (/o(m,0,/i)) 



8nL 



-mL|n| 



n 



(A.5) 



Next we consider the integral, /i(m. A, /i), which appears in the one- loop contribution to 
the nucleon mass, 



/i(m,A,/i) = fi' 
1 

~ 4 



{S ■ qf 



1 



(27r)" V ■ q — A + ie q^ — vrfi + ie 
Ai?i(m,^) + (A^-m2)/o(m,A,/i)) , 



where 



Ri{m, n) = jj'' J 



dJ'q 1 
(27r)" q^ — m? + ie 



167r2 



-m 



m 



r e + i-iog — 



(A.6) 



(A.7) 



This integral contributes to the pion-tadpole diagram in Fig. ^e). Subtracting the 7ii?T{e) 
divergence, one then has 



/i(m, A,/i) 



327r2 



(m^ - A^) \^/S? -m? + ie lo^ 

,2 \ 



'A - VA2 -m2 + ze^ 
A + VA2 -m? + ie 



m 



1 



- A log ^ + 2A 1 - -m^ A log — + -m^A 



327i' 



-F(m,A,/i) . 



Using the master relation, eq. ()A.3|) . one finds 



6l ( -Ri(m,/i) 



47r2L 



E 

n^O 



Ki{mL\n\) 



n 



(A.8) 



(A.9) 



Finally, we find [3 



(5l ( Ii{m,A,n)) 



167r2 



/5L|n| 



i^o(/3L|n|) 



(A.IO) 



Another useful integral is Jo{m, A, /i) = dIo{m, A, fi)/dA, 

(P^q 1 1 



Jo(m,A,/i) = 1^' J 



(27r)" (go - A + ie)2 g2 _ ^2 _^ 



log 



A 



/i2 VA2 - m2 + 



2e 



log 



'A - VA2 - m2 + ie' 



A + VA2 - m2 + ie 



(A.ll) 



The finite-volume corrections are 

iL 



5l(Jo("^, A,/i)) 



47r2 



d(3 



V/^2 + A2 - m2 



^ |n| iri(/5L|n|). (A.12) 
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The one-loop contributions to wavefunction renormalization, Fig.^f,g), and to the vertex 
diagrams, Fig. ^a,b), for the axial- vector current operator depend on 



Ji(m,A,/i) = /i' 
1 r 



(S-qf 



1 



4 



(27r)" {v ■ q — A + ieY q^ — fn'^ + 
i?i(m,/i) + 2A/o(m, A,/i) + (A^ - m^) Jo("^, A, /i) 



3 i 



4 167r2 



3 i 



m 



{m' - 2A') log — 



2Aa/A2 -m^ + ie Ioe 



'A - VA^ - m2 + ie' 
A + VA2 -m"^ + ie, 



J(m, A,/i) . 



4 167r2 

The finite-volume corrections may be written as 



h (Ji(m, A,/i)) 



167r2 



df3- 



E 



/3L|n| 



Ko{PL\n\) 



(A.13) 



(A.14) 



+ A2 - m^f ^0 

Finally, the vertex diagrams. Fig. ^c,d), for the axial- vector current operator depend on 

{S-qf 1 1 



iVi(m,A,/i) = /i^ 



(27r)" f-g — A + ie v ■ q + ie q^ — rn^ + ze 



— [ /i(m, A,/i) - /i(m, 0,/i) 



3 i 



4 167r2 



2 2 ^ 2m 



-|- -AV A^ — -|- ie loj 



'A - VA2 - m2 + ze' 
A + VA2 - m2 + ie^ 



3 i 



4 167r2 

The finite-volume corrections are simply 

1 



-|-^^ ( vrm — \/ A'^ — m? + it log 



K{m, A, yu) 



'A- VA2-m2 + 



le 



A + VA2 



(5L(Ari(m,A,/i)) = - [ 5i(/iKA,/i)) -5i(/iK0,/x)) ] 



(A.15) 



(A. 16) 



where one uses eq. ()A.10|) . 
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